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THE LINE CONTINUUM LUMINOSITY RATIO IN AGN; OR ON THE "BALDWIN EFFECT" 

R. MushotzKy and G.J. Ferland^ 

Laboratory for High Energy Astrophysics 
NASA/Goddard Space Flight Center 
Greenbelt, Maryland 20771 

ABSTRACT 

The luminosity dependence of the equivalent width of CIV In active 

galaxies, the "Baldwin" effect, Is shown to be a consequence of a luminosity 

dependent Ionization parameter, U. The available data, combined with a 

- 25 

"standard" photoionization model yield U a L ’ . This law also agrees with 

the lack of a "Baldwin" effect In Lya or other hydrogen lines. A fit to the 
available data gives a weak Indication that the mean covering factor decreases 
with Increasing luminosity, consistent with the Inference from X-ray 
observations. 

We discuss. In the appendices, the effects of continuum shape and density 
on various line ratios of Interest. 


^University of Kentucky 
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I. INTRODUCTION 

Observations of both high and low redshift quasars (Baldwin 1977; Baldwin 
et al. TD78) and Seyfert I galaxies (Wu, Boggess and Gull 1983) have shown 
that the equivalent width of the CIVM549 line In these objects Is lower for 
higher luminosity objects, the so called "Baldwin effect". Such a correlation 
Is Important because It offers an empirical method of determining the 
Intrinsic luminosities and distances of quasars (Baldwin et al. 1978) and may 
serve as a probe of the environment near the poorly understood central 
engine. Such an effect does not seem to be present In Lya although It may be 
present In Mgl 1X2798 (Baldwin et al. 1978). The reality of this "Baldwin" 
effect for radio quiet quasars has been questioned by Osmer (1980), even 
though It seems well established for radio loud objects. 

We wish to show In this paper that the Baldwin effect In radio loud 
active galaxies can be understood primarily as a systematic decrease In 
Ionization parameter with lu;/i1nos1ty and secondarily perhaps a decrease of 
covering factor with luminosity (although these correlations may not be the 
only "cause" of the "Baldwin effect"). Previous Interpretations have centered 
on the effects of Inbedded dust (see Shuder and MacAlpIne 1979). 

II. MODEL 

We use the photoionization code of Ferland (1981) In the manner described 
by Ferland and Mushotzky (1982). To calculate our models we have used solar 
abundances and an Ionizing Input spectrum of the form 

^ = A E exp -(E/Eq) + B E ph/cm^ sec eV 

ai =2.2 (Malkan and Sargent 1982) 

02 = 1.7 (Mushotzky 1982) 
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with A/B ~ 3 X 10^. This reprr^duces the mean spectrum of Seyfert I's quite 
well (Wu, Boggess and Gull 1983), and corresponds to a a of ~ 1.5, the mean 
value for X-ray detected optically selected AGN. We have chosen a cutoff 
value Eq a 200 eV so as not to exceed the observed soft X-ray flux (see 
Appendix A). The code uses methods and assumptions very similar to those 
described by Kwan and Krollk (1981) (see also Davidson and Netzer 1979). We 
win be primarily Interested In lines formed In the highly Ionized zone of 
the broad line clouds, where the physical conditions are less sensitive to the 
exact details of the model such as line trapping or X-ray effects which may 

strongly effect such lines as Ha, H3 or Mgll. As a consequence our 

theoretical results are basically similar to those found by Davidson (1977), 
Shuder and MacAlpIne (1979) and Welshelt, Shields and Tarter (1981). 

High luminosity quasars apparently have a Ic/wer X-ray to optical ratio by 
a factor of 50 (Reichert et al. 1982) than Seyfert galaxies such that 

A/B can be as large as 1.5 x 10^. However we find that the conclusions we 

reach about Lya, GUI] and CIV are relatively independent of the optical to 
X-ray ratio because these lines do not come from the "X-ray" heated deep zone 
(see Appendix A), however we do not examine the effect of possible variation 
In the UV index a, on line ratios. Because certain high Ionization lines like 
0 VI and NV and "low ionization" lines like Mgll and HB are more sensitive to 
the ratio of optical to X-ray luminosity, Lq/L^, one still should be careful 
In defining the continuum form. We have chosen models with constant density 
between log p = 9.5 and log p = 10.2 cm“^ and clouds with fixed column density 

The assumption of constant density rather than constant pressure seems as 
reasonable when it Is realized that the sound crossing time for typical 
clouds (6r ~ 10 cm, t 1 year) is of the order of the timescale for 
continuum variability. Therefore since the clouds are in virtually 
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instantaneous thermal and ionization equilibrium they cannot be in pressure 
equilibriuim. In our models as long as log is > 22.4 at/cm^, and log U < 
”1.0, the GUI], Lyot and CIV line fluxes are almost independent of how thick 
the cloud is at these densities and ionization parameters. In most of our 
model runs Cl 1 1] was saturated by a column density of log ~ 22.4, 
while Lya and CIV were saturated at lower N|^. 

Tlie model assumes that there is only one “slab" at a given distance from 
the central source, and spherical symmetry in the broad line region (see 
Davidson 1977). We use the same definition of ionization parameter U as in 
Ferland and Mushotzky (1982); U = Q(H)/47 t r^cng where Q(H) is the number of 
ionizing photons per second and thus U depends on the continuum shape where 
log U is typical of the order of -2 for active galactic nuclei (Davidson and 
Netzer 1979). 

III. RESULTS 

As is shown in Figure 1, at a fixed continuum luminosity, the 
CIII], Lya and CIV luminosities are a strong function of the density and 
ionization parameters. This has long been well known (see review by Davidson 
and Netzer 1979). Wu, Boggess and Gull (1983) point out that several objects 
with high CIII3/CIV ratios have high Lya/CIV and that several others which 
have low CIII3/CIV have low Lya/CIV. We wish to point out here that this 
correlation is very general. In Figure 2 we plot the ratio of CIII3/CIV 
vs Lya/CIV for all the AGN in the literature we know of. On this graph we 
also draw lines of constant log p but varying ionization parameter, log U, 
from our model calculations. Note that because of the strong thermostat 
effect of these cooling lines, their strength is not strongly sensitive to the 
carbon abundance assumed (cf. Shuder and Mac Alpine 1979). Virtually all of 
the points in this diagram lie between the lines log p = 10 and log p = 9.5. 
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Assuming tnat this represents "correct" physics we can restrict ourselves to 
the majority of points that cluster on the log p = 9,5 line. We can then 
assign an effective ionization parameter log for a given object 
appropriate for Ly«,cni3 ant' CIV, from the location of the object In Figure 
2. We have chosen those objects that cluster with ±10% of the theoretical 
line, consistent with an estimate of observational errors. The many fewer 
points that appear to duster around the other two lines do not allow a good 
estimate of the dependence of log on L to be determined. 

We show In Figure 3 for Seyfert galaxies and QSO's a plot of log 

vs log Li 45 q, the luminosity of the source at 1450 A . This latter value 
may be taken as Indicative of the total luminosity of the scarce as measured 
by Baldwin and thus Is appropriate for the Baldwin relation. It is quite 
clear from this plot that log Ugff Is inversely related to Log 11450 * This 
relation has a linear regression coefficient R = .90 and Is thus significant 
at greater than the 99.999% confidence level. The best fit power law gives 
log U « -.25 log L 1450 . The shall size of the exponent requires a large total 
range In L to see (which is why this effect was not noticed In samples 
restricted to high luminosity high redshift QSO's) but we feel It is quite 
significant. 

It Is Interesting to note that comparison of the mean of CIII3/CIV 
and Lya/CIV ratios for quasars and Seyferts separately, which has a much 
larger number of objects, also shows the same effect. The mean of Lya/CIV for 
high z QSO's from Osmer (1980) and Osmer and Smith (1981) Is -^2.5 while for 
Seyferts I's, of considerably lower luminosity, the ratio Is 1.5 (Wu, Boggess 
and Gull). For CIII3/CIV the mean ratio for QSO's Is .5 while for Seyferts It 
Is .2. Since the "mean" luminosity at 1450 A for the Osmer and Smith QSO's 
Is log - 31.5 while for the Wu et al. sample It Is log LJ 45 Q ~ 29.0, 
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this ;.(lso Indicates In our representation of the data a luminosity dependent 
Ionization parameter^ One must of course also remember that the CIII]/CIV 
ratio Is density sensitive and that some fraction of the observed variation 
could easily be due to variation In density from object to object {see 
Appendix B). However Lyot/CIV has only a very weak density dependence {for log 
U = -1.5, Lya/CIV ~ 9.0 < log p < 10.2 from ojir calculations) so that 

most of the Lyo/CIV variation cannot be due to density effects alone. 
Slmllarlly It Is well established that the Lya/CIV and CIII]/CIV ratios are 
Independent of the X-ray to optical ratio (Kwan and Krollk 1981) {see Appendix 
A) {whicn Is Itself also a function of absolute luminosity {Zamorani 1982)). 

We therefore argue, tentatively, that a major fraction of the observed 
variation In the Lya/CIV vs CIIIl/CIV plane Is due to variation In Ionization 
parameter and that this variation Is luminosity dependent. A similar 
suggestion for the origin of the Cl 1 13/CIV ratio being lower for Seyfert I's 
than quasars was made by Wu, Boggess and Gull {1983). 

IV. DISCUSSION 

A. Baldwin Effect 

Can this observed luminosity dependent Ionization parameter account for 
the Baldwin effect? Roughly speaking, the luminosity In a given line 1s LL^^g 
s (fl/4Tr) A|_fL(U) Iconic In this expression L^gp^. Is the luminosity In the 
continuum Integrated over the frequency range appropriate to Ionize the Ion 
emitting the given line, f|_{U) for several lines Is shown In Figure 1 and Al 
I s a constant for a given line. In the following discussion the covering 
fraction, n/4n, which Is the fraction of central continuum photons Intercepted 
by broad line clouds. Is set equal to unity. For CIV we have, for log U < 

-2.5 {Figure 1) that fL(U) ~ while for -2.5 < logU <- 1.5, 

fL(U) ~ Thus since U ~ if we associate L^q^^. with Lj^ 45 q one gets 
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kiV *-1450 ^eff ^ kiV ~ ^1450 ^ ^eff ^ 

-1.5. Using our calibration in Figure 3 of log 1^459 Log we see that 

log Ugff ~ -2.5 is achieved at log 13^459 ~ 31.0. This value of log 1^459 is 

appropriate for Baldwin's (Baldwin 1977, Baldwin et al. 1978) sample of high 

33+ 17 

redshift quasars. These authors find that iQjy ~ ^ 1459 ' compared to rur 
prediction of an index of .56 (the recent analysis of Kiang, Cheng and Zhou 
(1983) of >0smer's data gives gives L^jy a for radio quiet objects). 

We consider this a reasonable agreement. Of course this requires that 
be linearly related to 1 ^ 459 . A possible luminosity dependent covering 
fraction, such as has been observed in the X-ray band could make the agreement 
almost exact (see discussion in Section C.) More data points are necessary to 
define the log U vs Log L relation more clearly (Baldwin 1977, Baldwin et al. 
1978). 


Examination of Figure 6 of Wu, Boggess and Gull or Figure 6 of Wampler et 
al. (1983) shows that the Baldwin relation starts to fail badly at log 
Li 45 o ~ 29-29.5 (M 1559 = -22) or in our transformed units, cf. Figure 3, log 
^eff -1»75 more or less where fL(U) for CIV flattens (Figure 1). That is, 

we predict that the Baldwin effect should become much weaker for lower 
luminosity AGN. This occurs when fi_(U) a that is for Log L 1459 <30. 

even higher values of U and, thus lower values of L, f|_(U) is flat for CIV 
e.g. Lj.jy a U° and thus there should be ionization parameter “Baldwin 
effect" consistent with the recent results of Wampler et al . (1983). 

We thus feel that we can account for both the slope of the "Baldwin" 
relation and the place where it ceases to be significant from the L 1459 vs. U 
relation. The lack of a "Baldwin" effect in Lya, e.g. the fact that the 
luminosity in Lya is linearly related to the continuum flux. Is easily 
understood as well. Figure 1 shows that L[_y^ ~ U"*^®. Thus we predict 
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that ~ ^cont ^^f’s^stent with Osiner's (1980) results. 

B. Other Lines 

What are the other predictions of such a model? As Kwian and Krollk 
(1981) point out, the high Ionization lines such as OVI U035 and NV M240 are 
quite sensitive to log U variations. Our model would predict that low 
luminosity objects would have considerably stronger OVI and NV lines than high 
redshift quasars. Unfortunately there Is very little data on NV from Seyferts 
and none on OVI. 

If this Idea Is correct one might also see the "Baldwin effect" in an 
Individual variable object, or If It had the proper Ionization parameter 
perhaps even an "anti -Baldwin effect". An anti -Baldwin effect iii predicted 
for objects with a low luminosity since for large U, log U > -1.5, 
f(U) a U’*®. In fact, In NGC 4151, from the May 1978 to Oct 1978 lUE 
observations (Penston et al. 1981) show that the continuum at 2500 A decreased 
by 22% but that the flux In CIV did not change. On timescales much longer 
than the light travel time across the broad line region (.1-1 year) 
reverberation effects such as those discussed by Bahcall and Kozlovsky (1969) 
and Blandford and McKee (1982) are suppressed. Similarly In NGC 3783 (Barr, 
Willis and Wilson 1983) which has variability on a time scale of a month, 
longer than the light travel time across the broad llnjj region In the object, 
it has been reported that a Baldwin type relation holds. Fitting the Barr et 
al. data we find L^jy o» LJ^gQ. If In this object U a (we assume that 
the cloud position does not change), this gives using. Figure 1, a log 
Ueff ~ -1.75. Matching the variations In the Lyot/CIV vs CIII]/CIV plane by 
using Figure 2 for this object also give log Ugff -1.75. For a luminosity 
of log Lj^ 45 o « 27.0 this value of Log U Is In general agreement with the 
overall log L^^gQ vs. log relation (F1gure;3)* Thus In this one object 
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we get similar estimates for U from variations In line ratios^ from Its 
position In the continuum luminosity versus U plane and from the relationship 
between CIV and continuum luminosity. 

We therefore conclude that much of the available data are consistent with 
the "Baldwin effect" both In the ensemble of active galaxies and In N6C 3783 
being due primarily to luminosity related variations In the effective 
Ionization parameters. 

C. X-Ray Data and Covering Factor 

It has been claimed (Holt et al. 1980; Mushotzky 1982; Lawrence and Elvis 

1982) that the effective X-ray column density Is ^'elated Inversely to the 

X-ray luminosity If we assume that this Is a covering factor effect one 

might expect to see it In the Lyo and CIV line fluxes, that Is, since the 

higher luminosity objects have a lower covering factor, this would tend to 

depress the Lya and CIV line fluxes systematically with luminosity. 

Inspection of Figure 6 of Wu, Boggess and Gull shows that for log L < 29.0 

(where we predict that the Ionization parameter dependent Baldwin effect 

should be weak or not present) that Log LQjy l*® , at roughly 3a 

cont^ 

confidence, behaving as If the covering factor O was weakly dependent on L, 
e.g. n a L"’^. This small an effect would not be seen over a large dynamic 
range in L In the Lya or CIV line fluxes since It would be dominated by the- 
"Baldwin effect" at large L. This variation In covering factor Infered from 
the CIV line predicts that. In the X-ray energy range, one would see the mean 
covering factor change from a (postulated) value of 1 at ~ 42.5 (for NGC 
4151) to .6 at L^ ~ 43.5 which would totally account for the observed X-ray 
absorption effect seen by Mushotzky (1982). That Is, the fact that high 
luminosity active galaxies do not, typically, show X-ray absorption can be 
explained by a reduction In the mean covering fraction from 1 for low 
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luminosity objects like NGC 41S1 to ~ .7 for objects of a factor 10 higher 
X-ray luminosity. 

If we Include this luminosity dependent covering factor In our 
predictions, the slope of the best fit relation between CIV line flux end 
continuum luminosity changes from iQjy a to o *^cont* nearly 
perfect agreement with that measured by Baldwin and co-workers. We thus feel 
that all these results are compatible. It Is Interesting to note that If we 
extend this fit to the quasars of Log ~ 46 we would predict a covering 
factor of ,2, which Is similar to that Inferred for these objects by the Lyman 
continuum absorption arguments {Osmer 1980), 

D. Can the Baldwin Effect be Seen In the Balmer Lines? 

Inspection of Figure 3 shows that thi» f/ther hydrogen lines (Ha, H$) show 
similar weak dependences on log In particular, to compare with the data 

of Yee (1980) and Shuder (1981) we find that, roughly log L ~ -*4 log 
U, for Log U In the range -1.5 to -2.5 and log L ~ -*2 log U for log U < 
-2.25. Thus we would predict log L^^ -- 1.1 log Lgo^^ and log L ~ 1.05 
log LgQn^. The actual fit Yee gives for quasars Is log (L (Hp) ) ~ 1.00 log 
Lcont Shuder gives log (UHa) ~ 1.05 log L^gp^. The agreement between 
theory and observation again seems quite good. We have not taken possible 
systematic variations In covering factor (sec c) Into account here. 

V. CONCLUSIONS 

We have established that there exists a luminosity-ionization parameter 
correlation for active galactic nuclei (Quasars and Seyfert I's). This 
correlation co(,^,b1ned with photoionization model calculations predicts that the 
luminosity In CIV X1549 should be proportional to the continuum luminosity to 
the 1/2 power at high luminosity, e.g., the Baldwin effect. At low luminosity 
the effect should get weaker and disappear for low luminosity AGN's. 
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We can further show that such an effect should be present In Individual 
galaxies when their luminosity changes, as has been seen by Barr, Willis and 
Wilson 0983) for NGC 3783. 

In addition the calculations show that there should not exist a Baldwin 
effect in Lya, Ha, and H$ . Observational ly this has been shown by Yee 
(1980), Shuder (1981) and others. 

The lack of a Baldwin relation at low luminosity allows us to investigate 
the covering factor, o/4ir, vs luminosity relation in low luminosity AGN. We 
find that, n a L"*^. This weak anti -correlation is important enough to 
account for both the observed anti -correlation seen In the X-ray between 
luminosity and X-ray absorption as well as to predict the low covering 
fractions observed in high luminosity QSO’s. 

We thank the referee for his careful reading of the manuscript and his 


insightful comments. 
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FIGURE CAPTIONS 


Figure 1 - Relative line fluxes, FlIU), for Lycx, CIV CIII] (1909A), 

MgII{2798), Ha, H 0 and Hell (4686) vs- the Ionization parameter log U at a 
fixed density (log p = 10.0) and fixed column density log = 22.8 
at/cm^ for the continuum specified In the text 

Figure 2 - Log CIII]/CIV vs log Lya/CIV for a sample of active galaxies. 

The lines are theoretical line ratio's from photoionization models. The 
tic marks are at log = -1, -1.5, -2.0, -2.5, -2.75 for log p = 9.5, 

log Ugff = -1.5, -2.0, -2.5^ -2.75 and -3.0 for log p = 10.0 and at log 
Ugff * -1.5, -2.0, -2.5, -2.75 and -3.25 for log p 10.2. The data marked 
by an open symbol are from lUE (Wu, Boggess and Gull 1983), Green et al. 
1980; Ulrich and Bolsson 1983; and Malkan 1983). The solid dots are from 
Baldwin and Netzer (1978), the solid diamonds from Oke and Korycansky 
(1982). The open square Is the "mean value" for Osmers quasars and the 
star the mean value for Seyferts. 

Figure 3 - Log luminosity at 1450 A vs log Ug^f. The data are drawn from 
Figure 2 from the points that cluster along the log p = 9.5 line. 

Symbols have same meaning as In Figure 2. The line Is the best fit power 
law to the data of slope -.25. The Oke and Korycansky points have been 
converted to 1450 A continuua following Wampler et al. (1983). Not all 
points In Figure 2 have continuum L^^jq available In the literature and 
so do not appear In Figure 3. 


TABLE 1, LINE RATIOS FOR MODELS 1-5 


log U « -3.0 



Model 1 

2 

3 

4 

5 

Lya/CIV 

23.4 

22.8 

23.9 

24,8 

28.5 

CIII3/CIV 

1.86 

2.05 

2.16 

2.13 

2.0 

Ha/H3 

17.3 

13.1 

12.1 

15.5 

18.7 

MgII/H3 

21.6 

13.1 

11.1 

18.4 

2.47 


log U = -2.0 


Lya/CIV 

1.48 

1.33 

1.42 

1.27 

1.13 

CIII3/CIV 

.18 

.155 

,16 

.15 

.13 

Ha/H3 

3.9 

4.1 

5.52 

4.2 

3.5 

MgII/H3 

2.9 

1.76 

2.00 

3.03 

1.7 


log U = -1.00 


Lya/CIV 

1.04 

1.67 

1.92 

1.32 

.90 

CIII3/CIV 

.08 

.12 

.13 

.103 

.08 

Ha/H3 

1.41 

1.36 

1.33 

1.35 

1.68 

MgII/H3 

.58 

.316 

.22 

.415 

.32 
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APPENDIX A: CONTINUUM EFFECTS 

In this appendix we will consider the effects of varying the continuum 
form on the strong emission lines, Ly«, Ha, HP, CHI], CIV and Mgll. We 
caution the reader that since Ha, Hp and Mgll are generated mostly in the 
thick cool regions of the cloud whose physics is not yet totally understood 
that these results are primarily illustrative. We have constant density 
clouds of thickness log N^ = 22.8 and density log p =10.0. We pick 4 model 

-ai -02 „ 

spectra of the form F^ = A-v exp/-(E)/E(j) + Bv ph/cm'^ sec eV with = 
1.7, ai = 2.2 and vary Eq and A/B. 

Model 1 has A/B = 300, which corresponds to a^^^ =1.51 with E^ = «. 

Model 2 has E^ = 200 eV and A/B = 300. 

Model 3 has E^ = 200 eV A/B = 7200 which corresponds to a^^^ =2.05. 

Model 4 has E^, = 200 eV A/B = 2.3 or a^^^^ =1.28 

Model 5 has Eq = 200, A/B = .4, a^^^ = 1.00 

This range of a^^^ from 1.00 to 2.05 covers, roughly, the complete range 

of a as seen by Zamorani et al. (1981). A quick inspection of Table I shows 

OX 

that the line ratios are almost invariant with respect to the models. That 
is, within the limited range of continuum models we have chosen, the line 
ratios are almost invariant with respect to the exact form of the soft X-ray 
(.2-5 keV) spectrum or to the optical/X-ray ratio, l.rj largest changes 
between models is seen between models 3 and 4 the a^^ =2.05 and 1.28 models. 
Model 4 has a systematically higher ratio of Mgll/Hp. Our models are 
consistent with the observational results of Kriss (1982) on the Lya/CIV ratio 
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showing only a weak anti -correlation with 

At higher values of L^/Lq e.g. =1.0, Model 5, the log p = 10.0 
= 22.8 clouds become more Ionized all the way through and the nature 
models change since there 1s no longer a very cool "vasty deep". This 
strongly affects the Mgll, and H. lines. 


log 
of the 
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APPENDIX Bs DENSITY EFFECTS 

The referee has also suggested that we should clarify the use of CIII] as 
both a density and Ionization parameter diagnostic. In this section we 
consider the CIII]/CIV ratio for a constant continuum and column density but 
variable density and Ionization parameter. In Table 2 of Appendix II we show 
our results for a model with A/B = 7200 (e.g. low X-ray flux =2.05), and 
no cutoff and constant thickness log N^ = 22.8 for 3 different densities and 
Ionization parameters. 


Table 2 
CIII3/CIV Ratio 


log p 


log U 
-2.50 
- 2.00 
-1.50 


9.5 

.817 

.27 

.092 


10.0 

.420 

.129 

.060 


10.2 

.289 

.093 

.053 


As Is easily seen the results are degenerate. One can get the same values of 
CIIIl/CIV for all 3 densities and Ionization paramters. Thus one cannot use 
only the CIII]/CIV ratio to determine density or Ionization parameter. 


However as Figure 2 and Table 3 show, for a given continuum Lya/CIV is 
sensitive almost totally to Ionization parameter and not density. 
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TABLE 3: Lya/CIV Ratio as a Function of Density and log U 

log p 


log U 

9.5 

10.0 

10.2 

-2.50 

3.88 

3.28 

2.85 

-2.00 

1.26 

1.08 

.94 

-1.50 

.784 

.755 

.71 


Thus over fairly wide ranges In density and Ionization parameter one can 
estimate both log U and log p by using the two line ratios Lya/CIV and 


CIII3/CIV If the measurements have accuracies of better than 10%. 
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